Evolution of Rogue Waves in Interacting Wave Systems 
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Large amplitude water waves on deep water has long been known in the sea faring community, 
and the cause of great concern for, e.g., oil platform constructions. The concept of such freak 
waves is nowadays, thanks to satellite and radar measurements, well established within the sci- 
entific community. There are a number of important models and approaches for the theoretical 
description of such waves. By analyzing the scaling behavior of freak wave formation in a model 
of two interacting waves, described by two coupled nonlinear Schrodinger equations, we show that 
there are two different dynamical scaling behaviors above and below a critical angle Oc of the 
direction of the interacting waves below 6c all wave systems evolve and display statistics similar 
to a wave system of non-interacting waves. The results equally apply to other systems described 
by the nonlinear Schrodinger equations, and should be of interest when designing optical wave 
guides. 



I. INTRODUCTION 



The occurrence of rogue waves, i.e waves that are at 
least twice the size of the significant wave height of the 
surrounding waves, has long been a well known and much 
feared subject amon^professionals in the off-shore and 
sea faring business ([inil; 111) . The first truly scientific mea- 
surement of such a rogue wave was done at the Draup- 
ner oil platform in the North Sea off the coast of Nor- 
way on January 1, 1995 (4), where the so-called Draup- 
ner wave or New Years wave was observed. Since then, 
a large number of measurements, using different tech- 
niques, have been made in areas prone to such rogue 
waves, and the concept of rogue or freak waves is now 
scientifically well established [HI; 0; S S 0)- The find- 
ings show clearly that the rogue wave phenomenon is 
far more common than predicted by linear wave models. 
Both shallow and deep water wave behavior can to some 
degree be captured by simple nonlinear models, e.g. the 
Kortewe g-de Vries or nonlinear Schrodinger equations 
(NLSE) liO; [iil). For deep water waves, modified NLSE 
models, such as the Dysthe equation (|12l), have been pro- 
posed to more accurately characterize such waves. Simu- 
lations of the exact dynamical equations have shown the 
formation of freak waves in the nonlinear stage of the 
Stokes wave modulational instability (^13i) . Some of the 
main features of rogue waves are nevertheless still elu- 
sive in the simplified model descriptions. Hammack et 
al. (p^ studied the nonlinear interaction between large 
amplitude water waves in the laboratory and observed 
a dependence of the modulational instability on the an- 
gle between two large amplitude water waves. Recently, 
Onorato et al ([TBI) showed that interesting behavior arose 
if freak waves were modeled using nonlinearly interact- 
ing water waves. This model was later shown to give 
rise to wave pattern close to present rogue wave observa- 
tions (jl6|)- 



II. MODEL 

In this paper we present properties of a model of non- 
linearly interacting water waves and investigate to which 
degree the interaction promotes extreme waves, both in 
terms of the frequency at which waves of different height 
are occurring, and of the wave system's short and long 
term behavior. We will focus on a model for rogue waves 
based on two coupled NLSE (0; |l5|; ITBh . which permits 
large-scale simulations needed to bridge the gap to real- 
world wave statistics. The model system consists of two 
interacting waves A and 5, both traveling into a region 
of interaction. The wave system and its time evolution 
is illustrated in Fig. la. 

We analyze the statistical and spectral properties of 
a system of two nonlinearly interacting water waves by 
performing extensive numerical simulations. All mea- 
sured quantities are obtained from multiple realizations, 
each starting with different initial wave envelopes A 
and B by giving them a small amplitude, noise term. 
The evolution of the wave system is calculated from the 
slowly varying wave envelopes A and 5, which are re- 
lated to the actual wave surface r] diS r] = r]A -\- Vb 
with r]A = (1/2) A{r, t) exp {ikxX + ikyH — iu;t)-\- c.c, and 
Vb = {l/2)B{r,t) ex.p {ikxX — ikyy — iu;t)-\- c.c, where 
c.c. denotes the complex conjugate, kx and ky are the 
X and y components of the wave vector, and co is the 
wave frequency. The latter is related to the wave vector 
via the water wave dispersion relation uj = y/gR^ where 

g is the gravitational acceleration and = y^^^"^h^ is 

the magnitude of the wave vector. The time evolution 
of the wave envelopes are given by the coupled NLSEs, 
displayed in Fig. la. As we can see the last term dictates 
the coupling between the two waves with the coupling 
constant C = C{0) which ultimately is given from the an- 
gle 6. By setting ( = we obtain an uncoupled wave 
system. 

In short we have (0; M M): the group veloc- 
ity components Cx^y = ujkx,y/2hz'^^ the group dis- 
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FIG. 1 (a) An illustration of the wave system and the NLSE. Two waves A and B with the angle relative the dichotome 
are interacting and give collective wave dynamics in the blue region. The coupling constant ( — (^(0) dictates the level of 
interaction between the two waves A and B. (b) Time evolution of the variance for — 7i/8 and different values on 8. (c) 
a collapse of axy, the energy 8, the nonlinear coefficient ^ and time t for different values on and for the uncoupled (u-cpl) 
system. The right-hand vertical axis shows the kurtosis 72 of wave surface r]{x,y).{d) A collapse for = 7r/8,7r/7 and 7v/6 as 
in (c) but including the coupling coefficient as + 2^)^. (e) The different scaling behavior above and below Oc- t is the time 
to reach the second maxima of a^y for systems below Oc and the time to reach the only maximum of a^y for systems above Oc- 
The green plus signs (+) are, from left to right, for = 7r/8,7r/7 and n/G. The blue, crosses (x), are for = 7r/5 and 7r/4. All 
with 8 — 0.2. The gray line shows calculated times for different ^s using the scaling behavior in Fig. Ic and d. In all cases, the 
errors are smaller than symbols and the number of grid points Nx x Ny — 256 x 256 is used in the x and y direction. 



persion coefficients a = uo (2ky — k^) /8/^^, P = 
uo (2kl — ky) /8/^^ and 7 = —Zuokxky/An^^ the nonlinear 
coefficient ^ = ujk^ {2^ and finally, the coupling coefficient 

(kx — Ik). The vector components k^ and ky can be 
expressed from the angle between the wave vectors and 
the dichotome as kx — Kj cos d and ky = K sin 6. We have 
normalized the spatial coordinates and wave envelopes 
as r K~^r\ A i<i~^A'^ and B n~^B' ^ and time 
as t ^ Ui)~^t' to obtain dimensionless coupled NLSEs in 
terms of the dimensionless, primed variables, and where 
n and uo are eliminated. For convenience, we have omit- 
ted the primes below. 

We simulate the dynamics of the two interacting waves 
by using a pseudo-spectral routine — calculating the spa- 
tial derivatives in the spectral domain with periodic 
boundary conditions, and the fourth-order Runge-Kutta 
algorithm as time integrator. We use the grid sizes 
Ax = /S.y = 1 and the time step At = 0.5. The number 
of grid points are Nx x Ny = 256 x 256 in Fig. 1 and 
Nx X Ny = 128 x 128 in Fig. 2. As initial conditions, 
we set the A and B to constant amplitudes Aq and 
with Aq = Bq having values in the range 0.05-0.12. To 



seed the instability and generate multiple realizations for 
each parameter setting we add a small amplitude noise 
to A and B of the order Aq/IOO and Bq/IOO. 



III. RESULTS 

The total energy E = 

J {\A{x,y,t)\'^ + \B{x,y,t)\'^) dxdy is a con- 
served quantity. Locally the energy fluctuates, 
and the wave energy fluctuations of the system 
can be quantified by the variance cr^^(f,t) = 

^ / {\A{x,y,t)\^ ^\B{x,y,t)\^ - ef dxdy, where 
8 = E/L'^ is the average energy density. 

In Fig. lb we display the fluctuations (j'^y{8,t) for a 
system of two interacting waves A and B measured at 
different times t and with different energies £ using the 
angle = tt/S. We see that a'^y{£,t) initially increases 
with time for all £ and reaches a maximal value, then de- 
creases and eventually stabilizes at a value that is larger 
for larger values of £. Moreover, Fig. lb also indicates 
that a'^y{£,t) can be described by a function of only 
one variable. Specifically, we assume that the shape of 
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a'^y{£,t) can be expressed in terms of a dimensionless 
scaling function of only one variable, Fcr{X). We write 
this as aly{£,t) = £^ Fcr{{n^£)^{u:>tY). For dimensional 
reasons we then expect a = 2, as F^r is dimensionless, but 
leave this for the collapse plot. 

Fig. Ic shows a collapse plot of 9 = 7r/8 7r/4 
and an uncoupled (u-cpl) system (one collapsed curve 
for each case), with the values a = 2, 6 = 1 and 
c = 1. Hence, we can read out a scaling law of the 
form (jly{£, t) = £^ F^{^£t), since ^ = 00/^^/2 — the non- 
linear coefficient. We also see that as X oo, Fa{X) 
approaches the constant value 0.5, which yields the prop- 
erty a'^y{£,t) ~ 0.5£^ as t ^ oo. The scaling function 
tells us that the time evolution of the wave system display 
the same behaviour and is only scaled by the wave energy 
(and the nonlinear coefficient). The right-hand vertical 
axis in Fig Ic show the positive values of the excess kur- 
tosis of the wave surface 72 = M4(^)/cr(^)^ — 3, as this can 
be mapped onto our dimensionless scaling function since 
the variance of £ is proportional to the fourth moment of 
the wave surface and £ is proportional to the variance of 
the wave surface. This only holds when the fluctuations 
of the envelopes are large, and therefore we only display 
the the mapping for 72 > 0, which are obtained from 
processes with a wider-than-normal-distribution. We see 
that the coupled wave system reaches a maximum kurto- 
sis of 3 and that all wave systems finally approach zero, 
the kurtosis of a normal distribution. Similar results has 
been obtained for simulations of shallow water waves in 
crossing seas (11). 

Note that for larger angles = 7r/5,7r/4 and uncoupled 
wave systems the curves fall on top of each other, in 
contrast to 6 = 7r/8, tt/T and 7r/6 which are shifted. It 
was noted in Refs. (flBI : [l6l ) that there is a critical angle 
Oc = arctan(l/y2) ^ 0.615 rad ^ 35.3° at which the 
interaction changes character, making systems with < 
Oc more unstable due to overlapping instability regions 
and forming rogue waves at shorter time scales. This may 
explain the difference between = 7r/8, 7r/7, 7r/6 < Oc 
where the waves interact strongly, and = 7r/5,7r/4 > Oc 
where the waves behave as uncoupled (or single) wave 
systems. 

The space-independent harmonic solution is given 
as Q, Aeq = Boexp[-i{^B^ + ^C^Dt] and A^q = 
Bq ex.p[—i{£^BQ-\-2(AQ)t], where Aq and Bq are the initial, 
average, values of the wave envelopes. Since Aq = Bq^ 
the average energy density is £ = (Aq -i- Bq) = 2Aq^ and 
we see that A^q = Aq exp[i(^+2C)A§t] = Aq exp[(z/2)(^+ 
2Q£t]. So in the region of overlapping instabilities of the 
two waves, the time evolution might be scaled by {^-\-2() 
rather than just ^. We test this as before by performing 
a collapse for = 7r/8,7r/7 and 7r/6. In Fig. 2d we see 
that all three curves fall onto one and we can read out a 
new scaling as (jly{£, t) = f ^ F^[{^ + 2C)£t]. 

The coupled wave systems with = tt/S^tt/T and 7r/6 
have two local maxima, with the second maxima close to 
(^ -h 2()£t ~ 29. For systems of waves with no region of 
overlapping instability, when 35.3° < < 68.02°, there 



is only one peak of F^r, indicating that overlapping in- 
stability regions are needed to form the first peak. For 
= 7r/5, 7r/4 and uncoupled waves, the only maximum 
of axy is found at t ^ 19. Our assumption is that the 
second maxima for = 7r/8,7r/7 and 7r/6 and the only 
maximum for ^ = 7r/5 and 7r/4 are evolved by additional 
instabilities seen in both systems and represent the same 
point in their development, but is reached at different 
times du to a different scaling. 

By plotting the time it takes to reach the second max- 
ima for systems of overlapping instabilities and the only 
maximum for systems of no overlapping instabilities for 
different 6>s, the two curves should then if our assump- 
tion holds intersect at Oc and demonstrate that that the 
different scaling behaviour is separated by Oc. In Fig. 2e 
we plot t for the two different scaling behaviors and ad- 
ditionally show the measured values for the different ^s. 
We indeed see that the two curves intersect at the criti- 
cal angle Oc ^ 0.1967r rad as anticipated. This, together 
with Fig. Ic and Fig. Id tells us that 



C7ly{0,£,t) 



£^ F^ [^£t] ifO>Oc 
£^F^[{^^2C{0))£t] ifO<Oc 



(1) 



It was also noted in Refs. that counter- 

propagating waves and waves propagating at large an- 
gles, 68.02° < < 90° give rise to a new instability per- 
pendicular to their direction of propagation. However, 
this large angle regime is outside the scope of the present 
article and is left for future studies. 

Typical data from ocean waves give (T?), the wave fre- 
quency 0.09 Hz, uj = 0.56s~^, and hz = 0.033m~^. We set 
A = B = 0.1/ hz = 3m. In Fig. 2a we display the wave 
surface 7^, the power spectrum ^lAp, and the surface ele- 
vation probability density function P{r]) of an interacting 
wave system with = 7r/8. Note that only the positive 
part of the surface elevation probability density is dis- 
played, as the distribution is symmetric around zero. The 
snapshots are taken where the largest kurtosis is found 
(See Fig. Ic). This is exactly at the end of the exponential 
growth regime, with the wave envelopes A and B forming 
waves fronts extending in the y direction, traveling in the 
X direction. In Fig. 2b we look at the second peak of axy, 
(^ + 2Q£t ^ 30, of the coupled wave (upper), and com- 
pare it with the (Jxy-iociax.imdi^ ^£t ^ 20, of an uncoupled 
wave system (middle), also with = tt/S. This is where 
the two systems are showing the largest waves, despite 
the fact that the coupled system has a larger kurtosis 
in Fig 2a. The largest wave found in the coupled sys- 
tem is 21.24m, significantly larger than in the uncoupled 
case, 15.75m. The distribution P{r]) (bottom panel of 
Fig. 2b) is significantly wider in the coupled case, show- 
ing a tenfold, or higher, probability of generating waves 
of 15m or larger. The coupling re-distributes the wave 
energy and introduce additional wave modes in the power 
spectra as a result of the coupling to the 5- wave. 

Asymptotically (Fig. 2c), measured at (^ + 2()£t ^ 55 
for the coupled system and ^£t ~ 30 for the uncoupled 
system, the differences are diminished, which also was 
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FIG. 2 In (a) wave data from snapshots of a coupled system at the time point giving the maximum kurtosis 3.2). Top 
figure shows the surface elevation 77 (in meters) of the largest occurring wave from all realizations and their averaged power 
spectrum. Bottom panel shows the surface elevation probability density. In (b) the point in time giving the largest waves in 
the coupled (cpl) and the uncoupled (u-cpl) system; axy, + 2C)^t ^ 30 and ^St ^ 20, respectively. From top to bottom; the 
largest wave and a power spectrum of a coupled system, similarly fo the uncoupled system, and finally, their surface elevation 
probability densities, (c) the same as (b) but at a later time, showing the asymptotic properties of the surface elevations for 
the two systems, cr^y, + 2C)^t ^ 55 for the coupled system and ^8t ^ 30 for the uncoupled system. In (d) an illustration of 
where the snapshots are taken. We in all cases generate multiple realizations using 8 = 0.2, = tt/S and the number of grid 
points iV^ X iV^ = 128 x 128. 



indicated from the kurtosis in Fig Ic. Fig. 2d illustrates 
where the snapshots of Fig. 2a-c are taken. 



IV. DISCUSSION 

An increasing amount of data is currently being accu- 
mulated to improve rogue waves statistics and to map 
out areas prone to rogue waves . It is clear that such 
waves must be generated through a nonlinear mechanism. 
However, the treatment of coupled wave systems using a 
fully three-dimensional fluid dynamic setting is out of 
computational bounds. Therefore, although the model 
presented here does not take into account the full com- 
plexity of such wave interactions (12), it leaves the most 



essential nonlinear parts intact, and renders a statistical 
interpretation possible by performing large scale simula- 
tions. Another interesting aspect of the model presented 
here is its suitability for optical experiments. Recently, 
the occurrence of optical rogue waves in white light prop- 
agation has been observed Such direct optical ex- 
perimental studies of extreme waves can be used as a 
means for discriminating between different waves mod- 
els. Indeed, also here the propagation of light in opti- 
cal fibers gives an simplified picture of extreme water 
surface wave dynamics, but in essence retain the impor- 
tant properties of such rogue wave formation. Measuring 
rouge wave statistics of different wave guides could thus 
be performed to experimentally investigate the effect of 
the wave coupling more closely. The simplest example 
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would be to investigate the difference between reflective 
and absorbing boundaries in an optical flbre. 

In conclusion, we observe two different dynamical scal- 
ing behaviors separated by a critical angle 9c — the direc- 
tion of the wave systems relative to the dichotome. We 
also observe an increased probability of rogue wave oc- 
currence in a coupled two-wave system compared with 
an uncoupled wave system. The increase is twofold; the 
interaction increases the nonlinear focusing of wave en- 
ergy into larger waves, and additionally, the time needed 
to evolve large waves is decreased. The second part is 
particularly important when considering rogue waves at 
sea as wind speeds and directions vary significantly with 
time HO) • The local wind will align the waves in the 
wind direction, while turning winds are often associated 
with an inhomogeneous wind field making waves gener- 
ated in different regions being radiated into neighboring 
regions. 
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